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SPECTRAL ANALYSIS FOR THE EXCEPTIONAL A^-JACOBI EQUATION 


CONSTANZE LIAW, LANCE LITTLEJOHN, AND JESSICA STEWART 


Abstract. We provide the mathematical foundation for the Ym-Jacobi spectral theory. Namely, 
we present a self-adjoint operator associated to the differential expression with the exceptional Xm- 
Jacobi orthogonal polynomials as eigenfunctions. This proves that those polynomials are indeed 
eigenfunctions of the self-adjoint operator (rather than just formal eigenfunctions). Further, we 
prove the completeness of the exceptional Ym-Jacobi orthogonal polynomials (of degrees m,m-|- 
1, m -f 2,.. .) in the Lebesgue-Hilbert space with the appropriate weight. In particular, the self- 
adjoint operator has no other spectrum. 


1. Introduction 

The classical orthogonal polynomials of Laguerre, Jacobi, and Hermite are the foundational ex¬ 
amples of orthogonal polynomial theory. As shown by Routh in 1884 [20], but most often attributed 
to Bochner in 1929 [3], these three families of polynomials are, up to affine transformation of x, the 
only polynomial sequences satisfying the following two conditions: First, they contain an infinite 
sequence of polynomials {pn}^=o, where pn has degree n, such that for each n £ Nq, y = pn satishes 
a second order eigenvalue equation of the form 

p{x)y” + q{x)y' + r{x)y = Xy , 

where the polynomials p{x), q{x), and r(x) are determined by the corresponding differential ex¬ 
pression (Laguerre, Jacobi or Hermite). Second, each of the eigenpolynomials is orthogonal in a 
weighted space where the associated weight has finite moments. 

In recent years, there has been interest in the area of exceptional orthogonal polynomials, which 
presents a way to generalize Bochner’s classification theorem. The most striking difference be¬ 
tween classical orthogonal polynomials and their exceptional counterparts is that the exceptional 
sequences allow for gaps in the degrees of the polynomials. We denote an exceptional orthogo¬ 
nal polynomial sequence by using “Xm”, where the subscript m = \A\ denotes the 

number of gaps (or the codimension of the sequence). We require that the associated second order 
differential expression preserve the space spanned by the exceptional polynomials, but no space 
with smaller codimension. Consequently, the coefficients of the second order differential equation 
are not necessarily polynomial. Remarkably, despite removing any finite number of polynomials, 
the sequences remain complete in their associated space. 

Research in the area of exceptional orthogonal polynomials did not develop from a desire to gen¬ 
eralize Bochner’s theorem; rather, the exceptional polynomials were discovered in the context of 
quantum mechanics where researchers were looking for a new approach, outside of the classical Lie 
algebraic [11, 13, 18] setting, to solving spectral problems for second order linear differential oper¬ 
ators with polynomial eigenfunctions. In particular, they were discovered in [5, 7] while developing 
a direct approach [4] to exact or quasi-exact solvability for spectral problems. The first examples 
of these exceptional polynomials were introduced in 2009 by Gomez-Ullate, Kamran and Milson 
[5, 7], who completely characterized all Ai-polynomial sequences. Their result showed that the 

1 


2 


CONSTANZE LIAW, LANCE LITTLEJOHN, AND JESSICA STEWART 


only polynomial families of codimension one (in particular, having no solution of degree zero) satis¬ 
fying a second order eigenvalue problem are the Xi-Jacobi and Xi-Laguerre polynomials. Explicit 
examples of the X 2 families were given by Quesne [18, 19], who used the Darboux transformation 
and shape invariant potentials to find these new families. 

Higher-codimensional families, including the X^-Laguerre and Jacobi exceptional polyno¬ 
mial sequences, were first observed by Odake and Sasaki [17]. Further generalizations were observed 
regarding two distinct types of Jfm-Laguerre polynomials by Gomez-Ullate, Kamran and Milson 
[9, 10]. These JCm-Laguerre polynomial families do not contain polynomials of degree n € N for 
0 < n < m — 1. Furthermore, Liaw, Littlejohn, Milson, and Stewart [14] show the existence of 
a third type of Xm-Laguerre polynomials. The Type III Xm-Laguerre polynomial sequence omits 
polynomials of degree n € No for 1 < n < m. This new class of polynomials can be derived from the 
quasi-rational eigenfunctions of the classical Laguerre differential expression by Darboux transform 
as well as a gauge transformation of the Type I exceptional Xm-Laguerre expression. 

Following the discovery of exceptional polynomials, there has been a desire to study the properties 
of these polynomials more rigorously. The explanation for existence via the Darboux transformation 
of the higher-codimension X^-Jacobi and Xm-Laguerre polynomials and a remarkable observation 
regarding the completeness of the Xm-polynomial families was given by Gomez-Ullate, Kamran and 
Milson [6]. Gomez-Ullate, Marcellan, and Milson studied the interlacing properties of the zeros for 
both the exceptional Jacobi and exceptional Type I and Type II Laguerre polynomials along with 
their asymptotic behavior [10]. The properties of the Type III X^-Laguerre polynomials is studied 
[14]. 

Further studies regard the spectral analysis for the polynomials. The spectral analysis for the 
Xi-Jacobi polynomials (for m = 1, A = {0}) may be found in [15] along with an analysis of 
properties resulting from an extreme parameter choice, and for the Xi-Laguerre polynomials, the 
spectral analysis was completed in [2]. For all three types of the Xm-Laguerre polynomials, a 
complete spectral study is completed in [14]. 

In Section 2 we introduce the X^-Jacobi differential expression along with some properties. In 
Section 3 we then apply the Glazman-Krein-Naimark theory to obtain a self-adjoint operator asso¬ 
ciated to the differential expression, for which the corresponding domain contains the exceptional 
Xm-Jacobi orthogonal polynomials (see Theorem 3.2). Further, we show the completeness of the 
exceptional X^-Jacobi orthogonal polynomials (of degrees m, m -|- 1, m -|- 2,...) in the Lebesgue- 
Hilbert space with the appropriate weight (see Theorem 3.3). Summing up, we present the spectral 
analysis of the X^-Jacobi differential expression. 


2. Some Properties of the Exceptional Xm-JACOBi Expression 


We begin by summarizing some properties of the exceptional X^-Jacobi expression as described 
in [5, 10]. The parameters a and /3 are assumed to satisfy 


(2.1) a, /3>— 1, a-|-1 — m —^ {0,1,..., m — 1} , and sgn(Q:-|- 1 — m) = sgn/3 

in accordance with [10, Section 5.2], unless otherwise noted. 

The exceptional Xm-Jacobi polynomial of degree n > m, is given in terms of the classical 

'j 00 

Pu ' > by 

J k=o 




(- 1 )” 

a -|- 1 -|- n — m 


i(a + /3 + n-m + l)(x-l)Pi “ ^^(x)P,i"+-^^^(x) 

+ (a-m + l)P^-"-2’^i(x)P,i“+^’^"^^(x) . 
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The exceptional X^-Jacobi polynomials satisfy the second-order differential equation 

Ta,p,m[y\{x) = XnVix) 

for X G (—1,1), where the exceptional Jacobi differential expression is given by 
Ta,/3,m[y]i^) ■= (1 - x^)y'’{x) + - a -{P + a + 2)x-2{l - x^) (x))'^ ^ y'{x) 

(2.2) -I- ^(a — (3 — m + l)m - 2/3(1 - x) (x))^ ^ y{x) 

and Xn = — {n — m){l + a + /3 + n — m). 

In Lagrangian symmetric form, the exceptional X^-Jacobi differential expression (2.2) writes 
1 


Ta,/3,m[y]ix) — 


[Wa,l3,m{x){l - X^)y'{x))' 


bbQ!,/3,m(^) 

+Wa,0,m{x) (^rn{a - P - m + 1) - 2/3(1 - x) (log{Pj^'^~Xf^~^\x))^ ^ y{x) 
for X G (—1,1), where Wa,i3,m is the exceptional Xm-Jacobi weight function given by 

Wa,,^^rn{x) = ^ ^ ■ 

The restrictions on a and /3 ensure that Wa,y,m{x) has no singularities for x G [—1,1] and conse¬ 
quently, all moments are finite. 

The exceptional Xm-Jacobi polynomials \ Pjnn'^ r orthogonal with respect to the weight 

I ’ j n=m 

function Wa,i3,mix)- 

The eigenvalue equation Ta^i3^rn[y] = Xy does not have any polynomial solutions of degree n for 
0 < n < m — 1. Despite this fact, it is interesting that the exceptional X^-Jacobi polynomials 

1 f form a complete sequence in the Hilbert-Lebesgue space L^((—1,1); Wa,i3^m), defined 

by 

L^((-l, 1); Wa,i 3 ,m) := |/ : (-1,1) ^ C : / is measurable and J \ff Wa,0,m < ooj . 

3. Exceptional Xm-JACoei Spectral Analysis 

We follow the methods outlined in the classical texts of Akhiezer and Glazman [1], Hellwig [12], 
and Naimark [16]. 

The maximal domain associated with Ta^i3^m[-] in the Hilbert space L^((—1,1), Wa,i3,m) is: 

(3.1) A = {/ : (-1,1) ^ C|/, f G ACioc{-l, 1); /, T^Am[f] e 1), . 

The maximal domain A is the largest subspace of functions of L^((—1,1); Wa,0,m) for which 
maps into L^((—1,1); Wc^yg^m)- The associated maximal operator is 

Sl0,m ■■ C 1), W„,y3,m) ^ T'((-l, 1), bb«,y3,™) 

where p is defined by 

(3.2) := T^,0,m[f] 

:=A- 
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For /, G A, Green’s Formula may be written as 
(3.3) 


where [•, •](•) is the sesquilinear form defined by: 

[f, 9 ]{x) = Wa, 0 ,m{x){l - x^){f'{x)g{x) - f{x)g'{x)) 
f 1 — fl -l- 

^^ ^ ^ ^ , {nx)g{x) - f{x)g'{x)) {x G (- 1 , 1 )) 


p( OL l,/3 1) 

m 


(^)) 


(3.4) 
and where 

[f,g]{x) |^=ii:= [f,g]il)-[f,g]i-l). 

By the definition of A and the classical Holder’s inequality, notice that the limits 
[/>S'](- 1 ):= lim [/,S'] (a:) and [/,s] ( 1 ) := lim [/,s] (x) 

x—>■—1+ 

exist and are finite for each f,g G A. 

The adjoint of the maximal operator in L^((—1,1); Wa,/ 3 ,m) is the minimal operator, 

Slp,m ■■ C 1), VF,,;3,m) ^ 1), 

where is defined by 

f G ■■= {/ G A| [f,g] = 0 for all s G a} . 

We seek to find a self-adjoint extension ia -b^((~l) 1); l^a„a,m) generated by 

which has the exceptional Xm-Jacobi polynomials as eigenfunctions. In order to 

I ’ j n=m 

achieve this goal, we need to study the behavior of solutions at the singular endpoints x = — 1 and 
X = 1 so as to determine the deficiency indices and find the appropriate boundary conditions (if 
any). 

First, we obtain the deficiency indices via Frobenius Analysis. They depend on the values of the 
parameters a and jd. 

The endpoints x = — 1 and x = 1 are, in the sense of Frobenius, regular singular endpoints of 
the differential expression Ta^p^m[-] = 0. We first apply Frobenius analysis to the endpoint x = 1. 


By multiplying the exceptional X^-Jacobi expression by 


X — 1 
X + 1 


, we obtain 


X — 1 
X + 1 


with 


and 


) - Ky{x)) ^{x- lfy'\x) - {x - l)p{x)y'(x) + q{x)y(x) 

))'(x-l) 


, . /3 - a - (a +/3 + 2)x 

p(x) = -—- 2 

X + 1 




g(x) = - /? - m + l)m - 2/3 (^log(P^ " ^))) (x - 1 )^ . 

Evaluating the above equation at x = 1 yields the indicial equation 

0 = r(r — 1 ) — rp(l) + < 7 ( 1 ) = r(r + a). 
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Therefore, two linearly independent solutions to Ta^ii,m[v\ — ^nV = 0 behave asymptotically (near 
X = 1, e.g. on the interval (0,1)) like 

zi(x) = 1 and Z 2 {x) = {x — 1)““ 


near x = 1. 

For all allowable values of a and 13, 


Jo 


Wa,i3,mix) dx <oo; 


while 


\z 2 {x)\^ Wa,p,m{x) dx < OO 


only for —1 < a < 1. 

In a similar manner, multiplying the exceptional X^-Jacobi expression Ta^jS^miv] — J^a,p,my by 
(x + l)/(x — 1), results in an indicial equation 

r(r + /3) = 0; 

and two linearly independent solutions will behave (asymptotically) like 

yi(x) = 1 and y 2 {x) = (x + 1)“^ 

near x = —1. 

For all allowable values of a and (3, 

/■o 2 

/ \yi{x)\ Wa^i3^m{x) dx < OO] 


while 


f 2 

J \y2{x)\ Wa,/s^m{x) dx < OO 


only for —1 < (3 < 1. 

As a consequence, we have the following results. 


Theorem 3.1. Let Ta^p^miy] — ^a,y,m be the exceptional Xm-Jacobi differential expression (2.2) on 
the interval (—1,1). 

(1) is in the limit-point case at x = —1 for (3 > 1 and limit-circle for —1</3<1. 

(2) is in the limit-point case at x = 1 for a > 1 and limit-circle for —1 < a < 1. 

Corollary 3.1. The minimal operator has the following deficiency 

indices: 

(1) For a, 13 > 1, p has deficiency index (0,0). 

(2) For a > 1 and 0 < f3 < 1, has deficiency index (1,1). 

(3) Similarly, for P > 1 and 0 < a < 1, p has deficiency index (1, 1). 

(4) For —1 < a,/3 < 0, has deficiency index (2,2). 

Next we formulate the self-adjoint operators. 

Theorem 3.2. The self-adjoint operator Sci,y,m in L^((—1,1); Wa,0,m)^ generated by the exceptional 
Xm-Jacobi differential expression Ta^p^m is given by 

Sc.Am[f] = T^Am[f] 

f G 
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where 

(3.5) 




A 


if Oi>\ and (3 >\ 

{/e A 

lim (1 J- x)^~^^f'{x) = 0} 

x— >■—1 + 

if ot>\ and 0 < /3 < 1 

{/e A 

lim (1 — x)°~^^f'{x) = 0} 

X^l~ 

if0<a<l and f3 >1 

{/e Aj 

lim (1 — x)'^~^^f'{x) 

= lim (1 J- x)^~^^f'{x) = 0} 

x^—l+ 

for all other choices of parameters 
that are allowed by (2.1). 


Proof. If the parameters satisfy a, 13 > 1, then there is only one self adjoint extension (restriction) of 
the minimal operator (maximal operator 5*^^^ „j); that is, the maximal and minimal operator 

coincide and 

Suppose that 0 < a < 1 and f3 > 1, then there are infinitely many self-adjoint extensions of 
the minimal operator From Corollary 3.1, the deficiency index equals (1,1), which means 

that ^{3^ ^) is a subspace of codimension 2 in A. We will restrict the maximal domain A by 
imposing a suitable boundary condition which is invoked by the sesquilinear form [•,•](•) defined 
by (3.4). First note that h{x) = (1 — x)““ € A since 


Ta,j3,m[{l -x) “] = C>((1 - x) “) (near x = 1), 

which implies Ta^/3^rn[g] G because a < 1. Further, the constant function 

satisfies 1 € A and 


(3.6) 


[h, 1] \lzU= [h, 1] (1) = - 


a 


2 / 3+1 


p(-a-l,/3-l) 

m 


= a2^+^ / 0, 


where we used standard identities for the Jacobi polynomials and the Gamma function to hnd 

r(— a-|-m) T{f3 + m — a — l) T{—a + m) m!r(—a) 




ml r(/3 + m — a — 1) 


F(-a) 


m!r(—a) m!r(—a) 


= 1 . 


In particular, we obtain from equation (3.6) that the constant function 1 ^ 'D{S^^p^rn)- 

For 0 < /3 < 1 and a > 1, we can prove the corresponding statement in a similar manner. Lastly, 
for a,/3 < 1, we combine the above cases. □ 


Note that every polynomial, in particular the Am-Jacobi polynomials, will satisfy all of the 
boundary conditions given by (3.5). 

Next, we adapt ideas introduced in [8] and further developed in [14] (for the case of exceptional 
Laguerre orthogonal polynomial systems) to prove that the spectrum of the self-adjoint operators 
from Theorem 3.2 consists exactly of the eigenvalues corresponding to the exceptional Am-Jacobi 
polynomials (and nothing more). 

Theorem 3.3. The exceptional Xm-Jacobi polynomials form a complete set of eigen- 

L ’ j n=m 

functions of the self-adjoint operator Sa,/3,m in- L^((—1,1), hFa^^g^m)- Additionally, the spectrum 
(x{3a,p,m} of Sa,i3,m is pure discrete spectrum consisting of the simple eigenvalues 

= crp{Sa,i3,m) = {-{n - m){l + a + /j u - m) I n > m} . 
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Proof. The eigenvalue equations follow by the Darboux relations. 

It remains to prove the completeness of \ P^n^ r ™ I); Fix a,/3 in the 

L ’ J n=m ’ ’ 

allowed range, pick f € PL = L^((—1,1), Wa^p^m) and choose e > 0. 

Define the function 

fix) 


fix) ■■= 


Pr 


(-0-1,/3-1) 


From the relationship 




(x) 

WaAx) 


between the exceptional and the classical weight iWa^/s^m and Wa,/3) respectively), it easily follows 

WfWn = \\f\\L‘2{{-i,iy,Wa,p)- 

In particular, we have / € L^((—1,1); W^A- 
Next we apply Lemma 3.1 with the function 

and obtain the existence of p € P such that 

J-PA^-PP-^){x)pix) 

Let N be the degree of p. With this polynomial p we can compute 




> 


f - pA^-^’P-^\x)p{x) 


L2((-l,l);iy,,^) 


f - {pA"-^'^-^\x))\{x) 


H 


Our goal is to show that the approximant (^P^ ^ p(x) is contained in the (closure 

of the) vector space spanned by the exceptional Jacobi polynomials. To this end, we consider two 
(n + m + l)-dimensional vector spaces 

^n+ 2 m := [Pm’P ■- j = m, m + 1 ,..., n + 2 m| and 

J~n+2m {9 £ Pn+2m • (1 T Xi)Q (xj) + fiqixi) =01-, 

where we let Xi denote the m — 1 roots of the polynomial (x). The space LFn+ 2 m is 

motivated by the exceptional term in the exceptional Jacobi differential expression. Clearly, we 

have (^pA‘^~^'^~^\x)'^ p{x) G Pn+ 2 m- Since (imiFn+ 2 m = dim£’n+ 2 m we achieve our goal, if we 
can show that 

£n+2m C T” n+2m' 

Take Q G £n+ 2 m- Since £n+ 2 m is spanned by a basis of eigenvectors of the exceptional X^-Jacobi 
differential expression we have Ta,^A^n+ 2 m] C £n+ 2 m- It follows that 

To.,pAQ] ■= (1 - x^)Q" + (/? - a - (/3 + a + 2)x - 2(1 - x") (\og{PA"-^'^-^^)) Q' 

+ {a- (5 -m + l)m - 2/3(1 - x) (^log(P^“““F/3-i))^ q 
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is polynomial, and hence the exceptional term (that is, the only term with a denominator): 

-2(1 - x) ^ / [(1 + x)Q'{x) + l3Q{x)] 

Pm [x) 

is polynomial. Since the roots of the classical orthogonal are simple and 1 is not a root, we have 

(1 + Xi)Q'{xi) + PQ{xi) = 0 . 

We obtain Q G Pn+ 2 m as desired. □ 

Let P denote the set of all polynomials. 

Lemma 3.1. Given a function r] on [—1, 1] that satisfies 0 < c < \'i]{x)\ < C < oo for all x € [—1, 1]. 
Then the set {ri{x)p{x) : p G V} is dense in L^((—1,1); Woi,/ 3 ) for the classical range of parameters 
a,fi, and the classical Jacobi weight Wa,g- 

Proof. Fix a, fi in the classical parameter range; that is, a, fi > —1. Then, by the theory of classical 
orthogonal polynomials, the polynomials P are dense in Therefore, it 

suffices to show that 

P C cloS'^(? 7 'P). 

To show this, take p aP and fix e > 0. First observe that 

\\p/r]\\n < (l/c)||p||-H, 

so that p/rj ^TL. By taking q £P such that 

> C^Wp/v - qWu > \\{p/v - q)ri\\n = \\P “ mWu, 

the lemma is proved. □ 
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